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Abstract 

Phenological timing - i.e. the course of annually recurring development stages in 
nature - is of particular interest since it can be understood as a proxy for the climate 
at a specific region; moreover changes in the so called phenological phases can be 
a direct consequence of climate change. We analyze records of botanical phenology 
and study their fluctuations which we find to depend on the seasons. In contrast 
to previous studies, where typically trends in the phenology of individual species are 
estimated, we consider the ensemble of all available phases and propose a phenological 
index that characterizes the influence of climate on the multitude of botanical species. 
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1. Introduction 

Phenology is a well-known concept in ecology to describe the timing of cer- 
tain periodical development stages of species throughout the year fT]. Develop- 
mental stages, or phases (e.g. flowering, fruit ripening, leaf coloring, foliation), 
have been studied over many decades in Europe using defined plant species. 
This information is often used to develop phenological calendars and describe 
natural seasons 

Phenological phases are sensitive to temperature [1, 0] , and shifts of phases 
are often regarded as the first signs of a change in climate @, @, 01 • An average 
earlier onset of plant phases of 3.8 days per 1 °C increase over the last decades 
has been observed for Europe, with negative shifts for spring and summer phases 
and positive shifts for fall phases 

A well-known phenological record is the cherry blossoming in Kyoto, Japan, 
which has advanced by 7 days between 1971 and 2000 Q. It has been shown 
that the fiowering dates of closely related species in Japan have responded to 
chmate change in a similar way I10| . Nevertheless, early flowering plants deviate 
from this trend, showing larger advances due to warming than later flowering 
species, which could result in an ecological mismatch in the future. 
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The reaction of plants to climatic changes is non-linear and not uniform 



111 . Il2l Il3l | . It has been observed that while the correlation between air temper- 



ature and the onset of spring and summer plant phases is strong, the correlation 
becomes weaker for fall phases [3, 13, It is suggested that later in the year, 
other factors like water availabilitj/, nutrition and pollution gain in importance 
over the influence of temperature [8[. Moreover, the temporal and spatial vari- 
ability of phenological trends differs between plants and is strongest for spring 
phases [l5|. Differences in the phenolo gica l response to climate warming may 



also result from locally adapted species [16 1. 



Large uncertainties remain about the future development of phenological 
phases. Several studies concentrate on the influence of temperature and, by 
assuming a linear relation between rising air temperature and changes in the 



phenological cycle, extrapolate possible future changes [17|, |18|. Usually, this 



temperature sensitivity is analyzed by fin ding the best correlation for the pre- 
ceding months of an onset date, see e.g. figTli. 20l isj . While these previous 
studies concentrated on temperature responses of specific phases or groups of 
phases, no integrated approach assessing changes in the annual phenological 
cycle has been developed so far. We therefore propose a phenological index, 
which characterizes the annual phenological cycle by taking into account both 
the shift of spring phases and the shift of fall phases simultaneously. Following 
this approach, more general conclusions about climatic influences on phenol- 
ogy can be drawn since more data is used, implying better statistics, and an 
average prospect is obtained. The method is applied to the state North Rhine- 
Westphalia, Germany. 

The paper is organized as follows. In Sec. [2] we present our concept of a 
phenological index. The data this work is based on is described in Sec. [H 
The results of our analysis are given in Sec. S] in three subsections regarding 
fluctuations, the phenological index, and correlations between the index and 
temperature records. In the last Section we discuss the results and give an 
outlook. 



2. Method 

Phenological events are referred to as phases, since they take place on a 
specific day of the year and occur at a more or less regular pace. For the 
phase 0p,t, i.e. the day of the year when the phenological event p takes place in 
year t, and the average phase over all years, {(t>)p, we consider the phase anomaly 

fp,t = (f>p.t - {(f>)p , (1) 

where (■) denotes the average over time and (0) is defined by tan(</)) := i'^'" '^^ 



2l| . see [Appendix A[ Accordingly, ^pp^t is the anomaly record of the specific 
phenological event p. In the calculations, all phases (being originally a day of 
the calendar year) are transformed to the range < < 27r by — > where 
y = 365 or y = 366. 
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Our analysis is motivated by the following perception. In a year with ad- 
vantageous climatological conditions, spring phases occur earlier than expected, 
e.g. as observed in an early flowering of Forsythia. In addition, fall phases occur 
later than expected, e.g. as seen in a late leaf falling of Pedunculate Oak. In 
contrast, disadvantageous years lead to delayed spring phases and premature 
fall phases. In order to capture this effect, for a given year we study the phase 
anomaly (difference between actual phase and average phase over all years) ver- 
sus the corresponding average phase. In this representation, in advantageous 
years the anomalies of spring phases (located at the beginning of the year) will 
be negative and the anomalies of fall phases (at the end of the year) will be 
positive. We propose to use the statistical increase of the phase anomalies as 
a function of the average phase as a measure of how advantageous the climate 
of the corresponding year was for the ensemble of plants. Thus, separately for 
each year, we study the parameters of a corresponding linear regression model 
for ippt against {4>)p'- 

ip;^,=at{(b)p+/3t, (2) 

where at is the slope from the phase anomalies at year t and fit is the intercept. 
From the regression, we also obtain the root mean square deviations, a^p, which 
are given by the standard deviation around the regression. 

The linear fit, Eq. ([5]), to = (/> — ((/>) versus {(f>) provides the coefficients a 
and /3 (for simplicity we skip the indices). Together with Eq. ([1]) one obtains 
(eliminating ip) 

0= (0)(a + l) + /3, (3) 

which shows that a corresponds to a temporary change of frequency. Figure[TJa) 
illustrates that a positive slope a is related to a low frequency anomaly causing 
early phenological phases in spring and late phenological phases in fall (see also 
Fig. [2]). In the same way, (3 corresponds to a temporary phase shift, as illustrated 
in Fig. [TJb) - all phases appear before or after the average. 

We understand phenological processes to be triggered by an annual cycle 
that is a compound of all relevant climatological features. In general, such a 
cycle is unknown but we think of it as illustrated in Fig. [T] Once it passes 
a certain threshold, and its derivative has the right sign, such as positive for 
spring or negative for fall, a specific plant is activated and a phenological phase 
takes place, e.g. flowering in spring. 

In |Appcndix B| we show that the slope a is also associated with an increased 
or decreased cycle in such a way that the integral over an idealized annual 
cycle C is approximately proportional to a, as suggested by Fig. [ija). 

3. Data 

The study region. North Rhine- Westphalia (NRW), is the most populous 
state of Germany (« 18 million residents in 2008; 34, 070 km^ total area). Two 
types of landscapes can be found in NRW: the North German lowlands with an 
elevation just a few meters above sea level, and the North German low mountain 
range with elevations of up to 850 m. The lowlands comprise the Rhine- Ruhr 
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Figure 1: Idealized cycle of advantageous and disadvantageous phenological years as well as 
premature and delayed years. Illustration of the coefficients a and /3 from Eq. l(2]l. (a) Non- 
zero a, which is the slope of the linear fit to 0— (0) versus ((j>), constitutes a temporary change 
of the frequency. In the o > case this leads to earlier phases in spring and later phases in 
fall, (b) Non-zero /3, which is the intercept of the fit, constitutes a temporary phase shift, 
with overall earlier {j3 < 0) or later (/9 > 0) phases. Compare with Eq. ((Sjl. For illustration 
we use a shifted cosine as periodic function F. 

Area which is one of the largest metropolitan areas worldwide. These landscape 
features are also expressed by distinct types of climate. While in the lowlands 
the mean annual temperature is 10 °C with an annual mean precipitation of 
620 mm, in the mountainous regions the mean temperature is 5 °C and an annual 
mean precipitation of up to 1, 500 mm is common as measured between the years 
1961-1990 lH]. 

Onset dates of numerous phenological phases have been collected in Germany 
by the German Weather Services (DWD) for the past decades. Observations 
are carried out two to three times in a week, which determines the temporal 
accuracy of the dataset. Since 1951 data for over 159 phases has been observed 
at around 660 stations in NRW. Due to incomplete datasets, especially before 
1970, we have reduced the number of stations to those providing sufficient data 
for our purposes over the whole period from 1951 to 2006 (see [Appendix C| ). 
As agricultural phases are strongly influenced by agricultural practices as well 
as breeding, and show a weaker relation to temperature changes f2S] , they were 
not considered in this study. Thus, we analyzed time series data of 17 mete- 
orological and phenological stations in NRW for 75 phases for the period of 
1951-2006 (cf. Fig. in]). In order to investigate the effect of temperature on phe- 
nology, annual mean temperature records from the nearest climate station to 
each phenological station were further taken into account. Temperature records 
are based on observational data of the DWD and were partly interpolated . 
In [Appendix C| we list the phenological phases and climatological stations. 
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Figure 2: Examples of phase anomalies versus average phase in the case of Diilmen near 
Munster in North Rhine- Westphalia for the years (a) 1960-1964, (b) 1970, (c) 1980, and 
(d) 1990. The filled symbols represent the various phenological phases and the solid line is a 
linear fit through the data by least squares. 



4. Results 

4-1- Fluctuations 

In Figure [2] we show examples oi ip = (j) — (0) versus {(j)) , namely for the 
years 1960-1964, 1970, 1980, and 1990 at the station Diilmen. During winter, 
i.e. approx. (0) < 50 and (0) > 300, no phenological activity is recorded. While 
in 1960 [Fig. [UJa)] the phenological phases appear more or less as in average, in 
1961 spring phases occurred prematurely. In 1962 all phases were delayed and 
in 1963 and 1964 the spring phases only. 

In 1970 [Fig. [Hb)] spring phases occured late > 0) leading to a negative 
slope a. In 1980 [Fig. [2{c)] less phases were recorded but on the basis of the 
available data it seems to have been a rather normal year. In 1990 [Fig.[2fd)] the 
early phases appear prematurely {ip < 0) indicating good conditions in spring. 

Next we want to address how strongly the phenological phases fluctuate. In 



order to quantify these fluctuations we use the Rayleigh measure [2J, |21| 



(70 = v/(cos</-)2 + (sin</.)2 . (4) 

Here (■) is the average over time, separately for each phenological plant. If cj) 
is spread uniformly over the period then tr^ is close to because the averages 
of the trigonometric functions are very small. If 4> is fixed then a^h is 1. Thus, 
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<(|)> [days] <i|)> [days] 

Figure 3: Fluctuations according to Eq. (|4)l of phenological phases (a) in the case of Diilmen 
in North Rhine- Westphaha and (b) for all stations. In both cases the average is over all 
available years. Small correspond to large fluctuations. The solid line in the background 
of (b) illustrates the wavy pattern of the fluctuations (functional form oc sin((/>) cos((/<)). The 
numbers indicate some phenological phases that are spread relatively strongly. 1: Hazel, 
Corylus Avellana: flowering; 112: European Alder, Alnus Glutinosa: flowering; 114: Cornel 
Cherry, Cornus Mas: flowering; and 177: Wild Brier, Rosa Canina: fruit ripening. 



values close to indicate large fluctuations and values close to 1 small ones. We 
use this quantity since the standard deviation of an angle is not well defined. 
We would like to remark that is independent of the regressions in Fig. [2] 

Figure [3] shows the spreading versus the average phase {(j)) . The result 
for the example from Fig. [2] is depicted in Fig. ^a) . Two phenological phases 
have small values of cr^ and accordingly large spreading - which is due to mea- 
surement errors. Apart from that, most phases show > 0.95 and only early 
phases exhibit larger fluctuations (smaller cr^), compare with Fig. [2ja). 

In contrast, the cr^ obtained from all stations [Fig. [3][b)] look smoother 
and four phenological phases have rather small a^-values (large spreading). In 
general, a kind of wave pattern can be observed and is illustratively traced 
in Fig. mb): Spring phases exhibit larger fluctuations, early summer phases 
smaller ones, late summer phases again larger fluctuations, and fall phases again 
small ones. Calculating standard deviations, similar patterns have been found 
for 35 plant phases and 29 butterfly phases |15j] . Assuming a wave (see Fig. [1] 
and |Appendix Bj ) those phases with small fluctuations coincide with large slopes 
(or small negative slopes) of an idealized phenological cycle. 

Deviations from the curve could result from measurement inaccuracies, since 
for some phases (e.g. fruit ripening) the exact onset date is difficult to deter- 
mine. Another reason could be that those phenological phases with small fluctu- 
ation are triggered by a sharp change while the others with larger fluctuations 
typically occur in seasons when the trigger is not as sharp. In other words, 
small deviations of the phenological cycle barely influence those phases that in 
average occur when the idealized cycle has a large slope. Contrariwise, small 
deviations of the phenological cycle do affect phases that in average occur when 
the idealized cycle has a small slope. 
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Figure 4: Values obtained for the phase anomaly fit Eq. ^ as illustrated in Fig. [2] for the 
years 1951-2006 (Diilmen). (a) slope a (pheno-index), (b) intercept /3, (c) root mean square 
deviations from the fit cr^, and (d) number of phenological phases used for each year. The 
brownish area in the background of (a) corresponds to the standard error. 



4-2. Phenological index 

Systematically applying linear regressions to ip — cj) ~ {<j)) versus {(j)) of the 
example station Diilmen (solid lines in Fig. [5]) we obtain a set of quantities in 
Fig. m plotted against the corresponding year. In Fig. Ufa) and (b) we show 
the two fit coefficients, namely the slope, a, and the intercept, /3, respectively. 
As pointed out in Sec. [H the former indicates how advantageous a year is. We 
consider it as a phenological index (pheno-index). As can be seen, a fluctuates 
from year to year roughly in the range —0.2 < a < 0.2. 

Two additional quantities of interest are depicted in Fig. |4l(c) and (d) . The 
root mean square deviations from the fit in Fig. [21 a^p, capture how uniform the 
annual cycle is or how homogeneously the phenological phases respond to the 
climate variations. We find that except from an outlier in the year 2001 (due 
to a measurement error) the residuals are stable with values around or below 
10 days. Remarkably, the outlier does not seem to affect much the values of a 
and fi in 2001. Figure lUJd) shows the number of phenological phases considered 
in the specific years; this is the same number of points appearing in the panels 
of Fig. [21 Somehow - for the example station - up to 1970 constantly 43 phases 
were recorded per year. In 1971 only 11 values are considered but still a and /3 
seem to have reasonable values, supporting the robustness of the approach. 
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Figure 5: Correlations between pheno-index and annual mean temperature for the years 1951- 
2006 (Diilmen). (a) Comparison of annual mean temperature Tmean (open diamonds) and 
pheno-index a (filled circles), (b) Scatter plot of Tmean versus a for all years. For a better 
comparison, the records are normalized according to Eq. JSj and Eq. JgJ, respectively. 



4-3. Correlations between phenological phases and temperature records 

It is known that the temperature is an important chmatological element 
influencing the phenological timing, in particular at springtime [3|. Next we 
want to inspect, how the pheno-index (slope a) is related to the mean annual 
temperature T,„ean- Figure ^a) shows both, a as well as T,„oan measured at 
the closest chmatological station, nearby Billerbeck, which is situated less than 
20 km from Diilmen. In order to compare the two quantities, we have normalized 
both records to zero average and unit standard deviation: 

Tl-ncan = (rmcan " (Tmce^n)) / CTT (5) 

as well as 

a(") = {a~ {a))/a^ . (6) 

We find a fair agreement between the course of both quantities. However, 
from 1999 onwards the normalized temperature values are above the normalized 
pheno-index. By definition, due to continuity reasons, a cannot systematically 
deviate from zero. Thus, further research is required to reveal if this is a sys- 
tematic deviation or within the statistical fluctuations. 

In Figure [Hb) the quantities a and Tmoan are plotted against each other 
for each year. The correlation coefficient for this example is 0.60, which is a 
satisfying result considering the noisy data of Fig. [2] 

Figure [S] exhibits the area under investigation. North Rhine- Westphalia. We 
identify for each phenological station the closest temperature station and cal- 
culate the correlation value between the pheno-index and the associated annual 
mean temperature record. The locations of the stations and the resulting cor- 
relation values are depicted in the map. The correlations vary between 0.28 
and 0.64 and do not seem to show any systematic dependence on the position, 
indicating that micro-scale spatial chmatological conditions may be dominating 
the pheno-index. In addition, the correlation value of a and Tmean does not 
significantly depend on the amount of phenological data of the stations. 
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Figure 6: Area under investigation, North Rhine- Westphalia, and correlations between an- 
nual mean temperature and pheno-index. Since, in general, phenology (open diamonds) and 
temperature (open squares) are not measured at the same site, we identify the closest pairs 
and attribute the correlation value as a color-coded circle to the center between the two. The 
exemplary station Diilmcn illustrated in the previous Figures is indicated by the red solid line 
circle. 
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For several reasons we do not expect much larger correlation coefBcients for 
individual phenological stations. Certainly, the phenological phases are also 
influenced by other factors, in particular precipitation or sunshine duration. 
Precipitation has weak spatial correlations and could be responsible for micro- 
scale influences. Hence, including further information of this kind could improve 
the correlations. In addition, the spring phases are influenced by the beginning 
of the year and possibly by the end of the year before but certainly not by the 
end of the actual year. Thus, considering seasonal values or taking into account 
at least part of the previous year could also lead to stronger correlations. In 
addition, the fall phases are known to react in a less pronounced way to tem- 
perature 14 , 3, S| ^iid could contribute to noise leading to reduced correlation 
values. 



5. Discussion and Conclusions 

In summary, we analyze phenological records, characterize the fluctuations of 
the phases, and introduce a phenological index. We find that the spring and late 
summer phases exhibit the largest fluctuations while the early summer and fall 
phases exhibit the smallest fluctuations. This may be related to the derivative 
of the annual cycle, such as of temperature. By plotting (for each individual 
station) the phase anomaly against the average of each phase and applying a 
linear regression through all phases, we obtain a measure for how advantageous 
a specific year was for the ensemble of flora at the corresponding site. The 
slope represents a temporary change of frequency and the intercept a temporary 
phase shift. In addition, we show that the slope of such a fit is approximately 
proportional to the integral over the idealized annual phenological cycle. The 
advantage of our approach is that it characterizes the multitude of climatological 
factors influencing the entire phenological ensemble considered. It smoothes 
the volatile phenological data to a combined index which helps to detect and 
quantify impacts on life and life cycles. It can be applied even when records are 
incomplete, something which can cause problems when phenological phases are 
considered individually. Generally, a similar index can also be calculated for the 
phenology of fauna, but this is a task we leave for future studies. 

We compare the phcno-index with the annual mean temperature and find 
some agreement. The correlation value for various stations in the area of inves- 
tigation varies between 0.28 and 0.64 whereas from this study we do not find 
any systematic dependence of the correlations on the location. We conclude 
that additional factors influence the pheno-index. The response of phenological 
phases to changes in temperature has been found to vary spatially, for exam- 
ple being stronger in more northerly latitudes Q. Further, a more pronounced 
spring advance has been described for maritime western and central Europe 



compared to the continental east [17[. Temperature has been identified as the 
main driver of spring phases, followed by the photoperiod length. Other factors 
such as precipitation, nutrient availability and soil properties showed only minor 
effect in comparison to temperature [2^ . Further research is needed in order to 
figure out how our approach relates to these previous findings. 
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Climate change induced shifts in phenology could disrupt the chain between 
pollinator and plants [2^. Phenological plant phases are key stages of plant 
development - changes in their timing might influence other species. Thus, 
alterations of phenological phases could disrupt interactions among species, e.g. 



within food webs jl5| . see also |27| . Evidence from various species indicate an 



insufficient rate of phenological adaptation concerning food webs to changing 



climatic conditions [28|. 



It has been found that correlations between air temperature and fall phases 
are less pronounced [l3, 0, 0|. Thus, the application of the ideas of this work 
could lead to a better consideration of fall phases within phenological analysis. 
There is considerable interest on how phenology will be affected by climate 
change, particularly in the context of ecology or agriculture. The phenological 
index could be used as the basis of projections obtained from climate models to 
project the changes in phenology. 



Appendix A. Calculating the average phase {(p) 

We want to calculate the average phase, (</)), from a set of angles (t>i. In 
order to account for the cyclicality of the phase, we do not simply average the 
phases, but consider the Euler relation, e"^ = cos<^ + isin0, average the sine 
and the cosine separately, and use the relation tan 9 = . 

However, writing the inverse, (0) = arctan ^ ^ , is not precise since 
the arctan-function does not take into account the signs of numerator and de- 
nominator. Therefore, most programming languages provide the two argument 
function atan2, which properly calculates the angle, 

{4>) = atan2 ((sin0}, (cos0}) . (A.l) 



Appendix B. Relation between the pheno-index ex and the anomaly 
of the phenological cycle 

Although, the annual cycle of phenological advantage basically can have any 
periodic shape, we assume such a cycle has the form of a sine- wave (for simplicity 
we skip the indices p and t): 

C{(j)) = A sin(0 + \)+B, (B.l) 

where A is the amplitude, (f> — vt the phase (frequency ly), A the phase shift, 
and B an offset. We would like to remark that it would be more meaningful 
to express as a function of C, since the annual cycle triggers the phenological 
phases. Due to climate fluctuation, the cycle C deviates from the average annual 
cycle within a year as well as from year to year. Since such an idealized cycle is 
unknown, we study the phenological signals. 

Here we show that the slope a is associated with an increased or decreased 
cycle in such a way that the integral over C is approximately proportional to a, 
as suggested by Fig. [TJa). 
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The integral over one period of the average annual cycle vanishes when we 
drop the offset 



[Asin((0) + (A))]d(0) =0. 



(B.2) 



The quantities (j) and A are spread around {(j)) and (A), respectively (assuming 
A = const.), and are in general different from the averages. Using Eq. ([3]) in 
(jB.ip we express the integral over C as 



[^sin(((/.)(a + l)+/? + A))]d(0) 

r 

A 



ia- 


f 1 


A 




{a + 


1) 


A 




{a + 


1) 


(cos 


(- 


2A 





-cos(((/))(a + l)+/3 + A) 

/ J — TT 

[cos(7r(Q! + 1) + /3 + A) - cos(-7r(Q; + 1) + /3 + A)] 
[cos(7r(Q: + 1)) cos(/3 + A) - sin(7r(Q; + 1)) sin(;3 + A) 



(a + 1) 



(B.3) 



Since /3 + A«— 7r/2 (in order to match the seasons within the calendar year) 
the second term is sin(/3 + A) « —1. For a close to 0, the first term goes like 
—Tra and assuming a + 1 w 1 one obtains 



[v4sin(((/))(a + l) + /3 + A))]d((/.) 



2TrAa 



(B.4) 
(B.5) 



Therefore, we conclude that a, as the regression slope to 0— {(f)) versus (0), is a 
measure for the anomaly of the phenological cycle with respect to early spring 
phases and late fall phases and vice versa. However, the unchanged maximum in 
Fig.[lja) is not very realistic. Having in mind a temperature change, one would 
expect an overall vertical shift of the cycle and accordingly rather an anomaly 
reflected in the offset B. Nevertheless, since there are few phenological phases 
in winter (where no events are measured) and since the analysis is performed 
statistically, a can still be considered as a measure for the in- or decrease of the 
phenological cycle. 



Appendix C. Details on the phenological data 

The observational program of wild plants includes the following codes: 1- 
20, 64-74, 112-135, 175-178, 213-228, which are listed at http://www.dwd.de > 
Climate -|- Environment > Phenology > Observation programme > Wild plants. 
In order to have sufficient statistics, we filter the phenology data according to 
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Table C.l: Phenological and associated temperature stations 



phonological stations 


climate stations 


ID 


location 


PIK-ID 


location 


52334110 


Kevelaer, Klevc 


19183 


Weeze-Hees 


53331130 


Ziilpich, Euskirchen 


19006 


Euskirchen 


53341120 


Frechen, Erftkreis 


19107 


Pulheim-Brauweiler 


53371170 


Hennef, Rhein-Sieg-Kreis 


19113 


Hennef 


54110000 


Aachen (DWD), kreisfreie Stadt Aachen 


19004 


Aachen 


54352120 


Imgenbroich, Aachen 


19119 


Monschau 


55342130 


Billerbeck, Coesfeld 


19175 


Coesfeld 


55344110 


Diilmen, Coesfeld 


19177 


Billerbeck 


57331490 


Heidenoldendorf, Lippe 


20339 


Lage, Kr.Lippe-Hoerste 


57359110 


Exter, Herford 


15208 


Vlotho-Valdorf 


57391110 


Minden, Minden-Liibbecke 


15182 


Minden-Hahlen 


57412130 


Biihne, Hoxter 


20031 


Borgenb 


57421110 


Giitersloh, Giitersloh 


20022 


Guetersl 


58326170 


Warstein, Hochsauerlandkreis 


20259 


Warstein 


58397360 


Obernetphen, Soest 


20013 


Siegen 


58422410 


Wunderthausen, Siegen- Wittgenstein 


20291 


Berleburg, Bad- Wunderthausen 


59230000 


Witten-Stockum, Ennepe-Ruhr-Kreis 


19221 


Witten-Stockum 



the following criteria: (i) phenological phases with at least 3 entries for one 
station, (ii) years with at least 3 pairs of {4>), Lp, and (iii) stations with at least 
30 years of data, whereas the presence of the years 1951 and 2006 is required. 
The phenological and associated temperature stations are listed in Tab. IC.ll 
Daily temperature records have been averaged to annual resolution. Missing 
temperature data has been interpolated [23 |. 
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